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We present an analysis of inclusive f2° baryon production and decays in 230.5 fb _1 of data recorded 
with the BABAR detector, fij? baryons are reconstructed in four final states (Q~iv + , Q,~-k + tt°, 
f2~7r + 7T + 7r~, E~ K~ir + ir + ) and the corresponding ratios of branching fractions are measured. We 
also measure the momentum spectrum in the e + e~ center-of-mass frame. From the spectrum, we 
observe Q° production from B decays and in cc events, and extract the two rates of production. 
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PACS numbers: 13.30.Eg,14.20.Lq 



The (ess) is the heaviest weakly-decaying singly- 
charmed baryon. It has been observed independently in 
several decay modes by different experiments [l| and in a 
variety of production environments, including e + e" col- 
liders operating at the T(iS) resonance 00,13], photo- 
production [!, H, 0] , an d hyperon beams [J So far, B 
meson decays to have not been observed. Several dif- 
ferent mechanisms could contribute, principally weak de- 
cays of the following forms: b — > ccs (e.g., B~ — ► Cl^B c ); 
b — > cus (e.g., £>~ — > ); and 6 — * cwd (e.g., 

£?~ — * 7r _ ). Beyond the requirement to produce 
at least one ss pair during fragmentation, we would ex- 
pect these three types of decays to be further suppressed 
by the limited phase space, by |T« S | 2 , and by needing to 
produce a second ss pair, respectively. Theoretical pre- 
dictions for branching fractions of individual two-body 
contributions vary from 0(1(T 5 ) to 0(1CT 3 ) HEED]- 

In this letter, we present a study of the Q,® baryon, 
reconstructed in four decay modes: f2~7r+, f2~7r + 7r°, 
f2~7r + 7r + 7r~, and S~ K~ir + ir + [12]. We measure the ra- 
tios of branching fractions for these modes, normalizing 
to — * f2~7r + ). The previous most precise measure- 
ments of these ratios are from an analysis of approxi- 
mately 45 events from six f2° decay modes Q. We then 
measure the spectrum of the momentum in the e + e~ 
center-of-mass frame (p*) and observe significant produc- 
tion of baryons in the decays of B mesons. 

The data for this analysis were recorded with the 
BABAR detector at the Stanford Linear Accelerator Cen- 
ter PEP-II asymmetric-energy e + e~ collider. The de- 
tector is described in detail elsewhere . A total 
integrated luminosity of 230.5 fb _1 is used, of which 
208.9 fb _1 were collected at the T(45) resonance (cor- 
responding to 232 million BB pairs) and 21.6 fb _1 were 
collected 40 McV below the BB production threshold. 

Simulated events with the f2° decaying into the rele- 
vant final states are generated for the processes e + e~ — > 
cc -► il c X and e+e- -► T(4S) -► BB -► il c X, where 
X represents the rest of the event. The PYTHIA simula- 
tion package [3| is used for the cc fra gme ntation and 
for B decays to VL° C , and the GEANT4 Ig] package is 
used to simulate the detector response. To investigate 
possible background contributions, additional samples of 
generic Monte Carlo (MC) events are used, equivalent 
to 990 fb _1 for T(45) events (e+e~ -> T(4S > ) -> BB), 
plus 320fb _1 for cc continuum events (e + e~ — > cc) and 
340 fb _1 for light quark continuum events (e + e — > gg, 
q = u,d,s). 

The reconstruction of an f2° candidate begins by iden- 
tifying a proton, combining it with an oppositely charged 
track interpreted as a 7r _ , and fitting the tracks to a com- 
mon vertex to form a A candidate. The A is then com- 
bined with a negatively charged track interpreted as a 



K~ (tt~) and fit to a common vertex to form an (2~) 
candidate. For each intermediate hyperon (A, S~ , f2~) 
we require the invariant mass to be within 4.5 MeV/c 2 
of its nominal value (corresponding to approximately 4, 
3, and 3 times the detector resolution, respectively). We 
form 7T candidates from pairs of photons in the electro- 
magnetic calorimeter, requiring the energy of each pho- 
ton to be above 80 MeV and the combined energy to be 
above 200 MeV. We require the invariant mass of the ir° 
candidate, computed at the event primary vertex, to be 
in the range 120-150 MeV/c 2 . 

Each Q~ (2~) candidate that passes the requirements 
is then combined with one or three additional tracks that 
are identified as pions or kaons as appropriate. For the 
Sl~7r + 7T° final state, we also combine the hyperon and 7r + 
with a 7T°. The 0° candidate daughters are refit to a com- 
mon vertex with their masses constrained to the nominal 
values. From this fit we extract the decay vertices and 
associated uncertainties of the f2° and the intermediate 
hyperons, the four-momenta of the particles, and the fl® 
candidate mass. For each intermediate hyperon we re- 
quire a positive scalar product of the momentum vector 
in the laboratory frame and the displacement vector from 
its production vertex to its decay vertex. 

To further suppress the background, we compute the 
likelihood ratio C = UiPi^i)/ [Ui Pifa) + UiPf^i)] 
for each f2° candidate, where the index i refers to the like- 
lihood variables Xi, andpi(xi) are the probability density 
functions for signal (S) and background (B). For a given 
candidate C has a value between and 1. The likeli- 
hood variables Xi are the logarithm of the f2 _ or S~ decay 
length significance, which is defined as the distance be- 
tween the production and decay vertices divided by the 
uncertainty on that distance; the momentum of the 
or S~ in the e + e _ rest frame; the total momentum of the 
mesons recoiling against the f2 _ or S~ in the e + e~ rest 
frame; and, for the _ 7T + 7r° mode, the ir° momentum in 
the laboratory frame. These variables (particularly the 
decay length significance) cover the expected range effec- 
tively with a limited number of bins. The distributions of 
these variables for the signal hypothesis are derived from 
signal MC simulations, and for the background hypoth- 
esis from generic MC events in which contributions from 
real f2° are excluded. Separate distributions are used for 
each final state when measuring ratios of branching frac- 
tions, and for each momentum range when measuring the 
momentum spectrum. 

To measure the ratios of branching fractions, we re- 
quire that p* > 2.4 GeV/c in order to suppress combina- 
toric background. Since the kinematic limit for pro- 
duced in B decays at BABAR is p„ lax = 2.02 GeV/c, only 
ft® produced in the cc continuum are retained. We also 
require that the value of C for each candidate is greater 
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than a threshold Co, chosen to maximize the expected 
signal significance for a given final state based on simu- 
lated events. We perform an unbinned maximum likeli- 
hood fit to the mass distributions shown in Fig. [1] The 
signal lineshape is parameterized as the sum of two Gaus- 
sian functions with a common mean; the background is 
parameterized as a first-order polynomial. In the fits to 
the data, the signal yield is a free parameter; the widths 
and relative amplitudes of the two Gaussian functions 
are fixed to values determined from a fit to simulated 
signal events. The mean mass is also a free parameter, 
except for the S~ K~ tt + tt + final state where we fix it to 
the central value obtained in fig — > f2~7r + in order to en- 
sure proper fit convergence. The masses are found to be 
consistent with one another and with the current world 
average [l[ within uncertainties. 

The numbers of signal events are 177±16, 64±15, 25± 
8, and 45 ± 12 (statistical uncertainties only) for the final 
states 0~7r + , £l~7T + 7r°, f2~ ir + 'ir + ir~ , and S~ K~ir + ir + , 
respectively. These correspond to statistical significances 
of 18, 5.1, 4.2, and 4.3 standard deviations, respectively, 
where the significance is defined as v / 2A£ and A£ is the 
change in the logarithm of the likelihood between the fits 
with and without an Og signal component. The fitted 
yields are then corrected for efficiency, which is defined as 
the fraction of simulated signal events, generated in the 
appropriate p* range, that are reconstructed and pass all 
selection criteria. Including the loss of efficiency due to 
the A and f2~ branching fractions, we obtain efficiencies 
of (8.6 ± 0.6)%, (2.5 ± 0.3)%, (4.3 ± 0.4)%, and (4.7 ± 
0.5)% for the four final states, where the uncertainties 
include systematic effects and are partially correlated. 
The systematic uncertainties on, and corrections to, the 
ratios of branching fractions are listed in Table U and 
discussed further later. We measure the ratios to be 



B(n° c -> rrvr+Tr+TT-) 

g(ng -> z-K-iT + ir + ) 
B{n° c n-TT+) 



= 1.27 ± 0.31(stat) ± O.ll(syst), 
= 0.28 ± 0.09(stat) ± O.Ol(syst), 
= 0.46 ± 0.13(stat) ± 0.03(syst). 



We also measure the p* spectrum of fig in order to 
study the production rates in both cc and BB events. 
Only the fl~ir + final state is used. The same reconstruc- 
tion, optimization of selection criteria, and fitting pro- 
cedures described above are applied, except that no re- 
quirement on p* is made. Instead, the f2g candidates are 
divided into nine equal intervals of p* covering the range 
0.0-4.5 GeV/c. We again require £ > £q and compute 
the efficiency in each p* interval as before with simulated 
signal events. In the numerator of the efficiency we count 
events with measured p* in the appropriate interval, and 
in the denominator we count events with generated p* in 
that interval: this definition removes the slight broaden- 



u 




2.65 2.7 2.75 2.8 2.85 2.9 

Mass of Q.° c candidate (GeV/c 2 ) 



FIG. 1: The invariant mass spectra for candidates passing 
the selection criteria. The data are fit with a double Gaussian 
lineshape on a linear background. 



TABLE I: Systematic uncertainties on the ratios of branching 
fractions, where R 1 = B(Q,° C -» fr7r + 7r°)/B(fi° -» Sl~ir + ), 
R 2 = B(Q.° C -> n--ir + Tv + Tv-)/B(n° c -» fi-7r + ), and R 3 = 
B(Q,° -> E-K--k+-k+)/B{Q,° -> n~7r + ). 



Effect 



Ri Ri 



R 3 



Finite MC sample size 0.7% 0.7% 0.8% 

Intermediate resonances in Q° decay 1.3% 2.6% 3.7% 

Signal lineshape 1.0% 1.0% 1.0% 

Dependence on the fit procedure 1.5% 1.5% 1.5% 

Hyperon branching fractions 1.0% 

Particle identification efficiency 1.0% 1.0% 1.0% 

Tracking efficiency 0.0% 2.8% a 2.8%° 

p* spectrum mismodeling 1.5% 0.6% 3.5% 

7T° fitting and efficiency 7.8% b 



Total systematic uncertainty 



8.3% 4.4% 6.3% 



°A relative correction of +0.5% applies to R2 and R3. 
b A relative correction of +1.1% applies to R±. 
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FIG. 2: The background-subtracted and efficiency-corrected 
fic p* spectrum. The black points represent the data, with 
vertical error bars giving the sum in quadrature of statisti- 
cal and uncorrelated systematic uncertainties. The solid his- 
togram shows the Bowler fragmentation function, binned and 
fit to the data for p* > 2 GeV/c (vertical, dashed line). 



ing effect of the detector momentum resolution. We also 
take into account a small difference in efficiency between 
cc and BB events. The efficiency-corrected yield in each 
p* interval is shown in Fig. [2] 

The systematic uncertainties are divided into two cat- 
egories: normalization effects, which are treated as fully 
correlated between all p* intervals, and shape effects, 
which are treated as uncorrelated between different p* 
intervals. The normalization uncertainties are due to the 
mass resolution, which is determined from the MC and 
checked with studies of the control modes S° — > S~7r + 
and S+ — > S~7r + 7r+ (2.4%); the A and Vt~ branching 
fractions []J (1.3%); and the tracking efficiency, which is 
corrected for data/MC discrepancies with control sam- 
ples of t decays (5.9%). The shape uncertainties are due 
to the limited size of MC samples (< 1%); dependence 
on the fit procedure (1.5%); modeling of the p* spec- 
trum, which can affect the weighted average efficiency 
within a p* bin (0-6%); the signal lineshape parameter- 
ization (1.0%); and the particle identification efficiency 
(2.0%). When fitting fragmentation functions (see be- 
low), we consider only the statistical and shape uncer- 
tainties, added in quadrature. When quoting total yields 
and rates, we include the normalization uncertainties, 
along with a relative correction of +1.0% due to a known 
data/MC discrepancy in tracking efficiency. 

The double-peak structure seen in the p* spectrum is 
due to two production mechanisms: the peak at lower 
p* is due to 0° production in B meson decays and the 
peak at higher p* is due to 0° production from the cc 
continuum. This is consistent with the pattern observed 
in A+ and S° spectra measured for e + e~ annihilation at 
y/s = 10.6 GeV [H, [13, El] ■ We nt the P* spectrum with 



the Bowler fragmentation function [lj| for p* > 2 GeV/c. 
We then extract the continuum yield as the sum of the 
data points above 2 GeV/c plus the integral of the extrap- 
olated function below 2 GeV/c. Similarly, the yield from 
B decays is the sum of the data points below 2 GeV/c 
minus the integral of the extrapolated function below 
2 GeV/c. Note that we do not fit a fragmentation func- 
tion to the data below 2 GeV/c. We obtain yields of 
2583 ±289 and 2426 ±414 for ft° produced in the contin- 
uum and in B decays, respectively, where the uncertainty 
includes all statistical and experimental effects. An addi- 
tional model uncertainty arises from the extrapolation of 
the continuum tail for p* < 2 GeV/c. To estimate this, 
we repeat the p* spectrum fit and yield measurement 
with other fragmentation functions: Collins and Spiller 
(CS) [2(| , two versions of the phenomenological model of 
Kartvelishvili et al. (KLP-M and KLP-B) 0, [H and 
the Peterson model (23[. The CS and KLP-M fits are 
inconsistent with the data for p* > 2 GeV/c. The RMS 
of the yields from the three other fits is 240 events and 
is taken as the model uncertainty for the B and contin- 
uum yields. Dividing the yield in B decays by the 
total number of B mesons in the data sample, we obtain 
the branching fraction product B(B — ► f2°A),B(f2° — > 
Q~n + ) = [5.2 ±0.9 (exp)±0.5 (model)] x 10~ 6 , where X 
represents the rest of the B meson decay products. Di- 
viding the yield from the continuum by the integrated 
luminosity and correcting for the small variation in cross- 
section with -y/s, we obtain the cross-section product at 
y/a = 10.58 GeV: a{e+ e - -> Q° c X)B(n° c -> 0-tt+) = 
[11.2 ± 1.3 (exp) ± 1.0 (model)] fb, where X represents 
the rest of the event. As a cross-check, we also make 
model-independent estimates of the yields from the con- 
tinuum and from B decays by subtracting the data below 
the T(45) threshold. Within large uncertainties, these 
are consistent with the yields measured above. 

It is thus clear that decays of B mesons to il" occur 
at a significant rate. Assuming the absolute branching 
fraction B(n° c -► fJ-7T+) - 1%, we conclude that B(B -> 
f2°X) ~ few x 10 -4 . This is substantially lower than the 
inclusive B meson branching fractions to the charmed 

baryons A+ and H c , which are ~ few x lo- 2 [HEIElI. 

One possible explanation for this is that both A+ and E c 
can be produced in a b — > ccs transition without creating 
an ss pair from the vacuum, whereas at least one ss pair 
must be created for production. It is also possible that 
phase space suppression in B decays to baryons becomes 
significant when very close to threshold. 

In conclusion, we have studied the baryon at BABAR 
through four hadronic decay modes, using 230.5 fb" 1 of 
data. We measure the ratios of branching fractions for 
four modes, significantly improving upon the previous 
values 0. We have also measured the p* spectrum and 
found comparable production rates of fl® baryons from 
the continuum and from B meson decays. The inclusive 
B branching fraction to f2° is found to be substantially 
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lower than those to an d A+ baryons, assuming the 
relevant baryon weak decay branching fractions are of 
the same order of magnitude. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput- 
ing organizations that support BABAR. The collaborat- 
ing institutions wish to thank SLAC for its support and 
kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
INFN (Italy), FOM (The Netherlands), NFR (Norway), 
MIST (Russia), MEC (Spain), and PPARC (United 
Kingdom). Individuals have received support from the 
Marie Curie EIF (European Union) and the A. P. Sloan 
Foundation. 



* Deceased 

t Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

* Also with Universita della Basilicata, Potenza, Italy 

§ Also with IPPP, Physics Department, Durham Univer- 
sity, Durham DH1 3LE, United Kingdom 
[1] Particle Data Group, W.-M. Yao et al., J. Phys. G 33, 1 
(2006). 

[2] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 

288, 367 (1992). 
[3] CLEO Collaboration, D. Cronin-Hennessy et al, Phys. 

Rev. Lett. 86, 3730 (2001). 
[4] CLEO Collaboration, R. Ammar et al, Phys. Rev. Lett. 

89, 171803 (2002). 



[5] E687 Collaboration, P. L. Frabetti et al., Phys. Lett. B 

300, 190 (1993). 
[6] E687 Collaboration, P. L. Frabetti et al, Phys. Lett. B 

338, 106 (1994). 
[7] FOCUS Collaboration, J. M. Link et al, Phys. Lett. B 

561, 41 (2003). 
[8] WA89 Collaboration, M. I. Adamovich et al, Phys. Lett. 

B 358, 151 (1995). 
[9] P. Ball and H. G. Dosch, Z. Phys. C 51, 445 (1991). 
[10] S. M. Sheikholeslami and M. P. Khanna, Phys. Rev. D 

44, 770 (1991). 

[11] V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B345, 
137 (1990). 

[12] The use of charge conjugate modes is implied throughout. 
[13] BABAR Collaboration, B. Aubert et al, Nucl. Instr. 

Methods Phys. Res., Sect. A 479, 1 (2002). 
[14] T. Sjostrand et al, Comput. Phys. Commun. 135, 238 

(2001). 

[15] S. Agostinelli et al, Nucl. Instr. Methods Phys. Res., 

Sect. A 506, 250 (2003). 
[16] Belle Collaboration, R. Seuster et al, Phys. Rev. D 73, 

032002 (2006). 

[17] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 75, 
012003 (2007). 

[18] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

95, 142003 (2005). 
[19] M.G. Bowler, Z. Phys. Cll, 169 (1981). 
[20] P.D.B. Collins and T.P. Spiller, J. Phys. G 11, 1289 

(1985). 

[21] V.G. Kartvelishvili, A.K. Likhoded, and V.A. Petrov, 

Phys. Lett. B 78, 615 (1978). 
[22] V.G. Kartvelishvili and A.K. Likhoded, Sov. J. Nucl. 

Phys. 29, 390 (1979). 
[23] C. Peterson, D. Schlatter, I. Schmitt, and P.M. Zerwas, 

Phys. Rev. D 27, 105 (1983). 



